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Abstract 



The Teukolsky equation has long been known to lead to divergent integrals 
when it is used to calculate the gravitational radiation emitted when a test 
mass falls into a black hole from infinity. Two methods have been used in the 
past to remove those divergent integrals. In the first, integrations by parts 
are carried out, and the infinite boundary terms are simply discarded. In the 
second, the Teukolsky equation is transformed into another equation which 
does not lead to divergent integrals. The purpose of this paper is to show that 
there is nothing intrinsically wrong with the Teukolsky equation when dealing 
with non-compact source terms, and that the divergent integrals result simply 
from an incorrect choice of Green's function. In this paper, regularization of 
the Teukolsky equation is carried out in an entirely natural way which does 
not involve modifying the equation. 
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I. INTRODUCTION 



In 1971, Davis, Ruffini, Press, and Price M computed for the first time the amount of 
energy carried off by gravitational waves when a test mass, released from rest at infinity, 
falls radially into a Schwarzschild black hole. By numerically integrating the recently derived 
Zerilli equation p[ for black-hole perturbations, they found an amount AE = 0.0104/i 2 /M, 
where \x is the mass of the infalling particle, and M the mass of the black hole. (Units are 
such that G = c = 1.) 

This calculation was later generalized to many different situations. In 1973, Ruffini 
H considered a particle released from infinity with nonvanishing initial velocity. In 1979, 
Detweiler and Szedenits |3J] examined infall trajectories with nonzero angular momentum. 
The first calculation involving a Kerr black hole was carried out in 1982 by Sasaki and 
Nakamura ||, who considered a particle falling along the symmetry axis. In 1983, Kojima 
and Nakamura || examined the case of a particle moving radially within the equatorial plane. 
In 1984, the same authors [[!]] studied infall trajectories with nonzero angular momentum. 
The infall of spinning test masses was considered for the first time very recently, by Mino, 
Shibata, and Tanaka ||. The gravitational waves emitted by a particle in bound motion 
around a black hole have also been extensively studied |9|-p6f , following the pioneering work 
of Detweiler in 1978 0. 

Davis et. al []J obtained their classic result by integrating the Zerilli equation [Q, 
which describes in a compact form the (even-parity, or polar) metric perturbations of the 
Schwarzschild black hole. On the other hand, Detweiler and Szedenits |4]] worked with the 
Teukolsky equation f28fl , which describes in a compact form the curvature perturbations of 
the Kerr black hole. (In the limiting case of a nonrotating black hole, and in the absence of a 
source for the perturbations, the Teukolsky equation reduces to the Bardeen-Press equation 
p9| , which was derived earlier.) 

However, when integrating the Teukolsky equation, Detweiler and Szedenits encountered 
divergent integrals, which they regularized by integrating by parts and then discarding the 
infinite boundary terms. (The Zerilli formalism does not lead to divergent integrals.) This 
procedure was also adopted by Simone, Poisson, and Will [^DJ who, in 1995, reproduced the 
Davis et. al result using the Teukolsky equation, for the purpose of testing the accuracy of 
post-Newtonian methods. In neither of these papers was a justification given for discard- 
ing the infinite boundary terms, and confidence in the procedure came entirely from the 
comparison with Davis et. al. 

The first attempt to remove the divergent integrals from the Teukolsky formalism came 
in 1981 from Tashiro and Ezawa [31], wno employed the clever trick of subtracting from the 
dependent variable a quantity constructed from the source. The resulting equation for the 
new dependent variable leads to well defined integrals. Also in 1981, Sasaki and Nakamura 
dealt with the divergent integrals in a different way. In the restricted context of a 



Schwarzschild black hole, these authors rewrote the Teukolsky equation in a Regge- Wheeler 
form, with a source term constructed from the source of the Teukolsky equation. (The 
Regge- Wheeler equation describes the odd-parity, or axial, metric perturbations of the 
Schwarzschild black hole.) The Sasaki-Nakamura equation, which was later generalized to 
the case of a Kerr black hole |5[], does not lead to divergent integrals. Subsequent studies of 
infalling particles were carried out by integrating this equation. 
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Those methods of regularization give the impression that when dealing with unbounded 
particle trajectories, the Teukolsky equation necessarily leads to divergent integrals, and is 
therefore not well posed. The purpose of this paper is to show that this impression is based 
on a misconception. Indeed, I wish to show that the Teukolsky equation can be regularized 
in an entirely natural way which does not involve modifying the equation. 

The basic issues are easily summarized. 

For perturbations of a Schwarzschild black hole, to which I specialize in this paper, 
and after separation of the variables (the usual Schwarzschild coordinates are used), the 
(inhomogeneous) Teukolsky equation takes the form 

dd +p2u \ R = p2T 
ar ar J 

Here, R[f) is the radial function corresponding to a perturbation of frequency u and 
spherical-harmonic indices £ and m, p(r) = (r 2 — 2Mr)~ 1 , U(r) is the effective potential, 
whose explicit expression can be found in Eq. ( [2.15| ), and T{r) is the source term, which is 
constructed from the energy-momentum tensor of the infalling mass. 

The inhomogeneous Teukolsky equation is solved with the physical requirement that 
gravitational waves must be purely ingoing at the black-hole horizon, and purely outgoing 
at infinity; this is equivalent to a no-incoming-radiation initial condition. Mathematically, 
this translates into two statements. First, that R(r) oc R H (r) near r = 2M. Here, R H (r) 
is a solution to the homogeneous equation, normalized such that R H (r — > 2M) ~ (1 — 
2M/r) 2 exp(— iur*), where r* = r + 2Mln(r/2M — 1). Second, that R(r) oc R°°(r) near 
r = oo, where R°°(r) is also a solution to the homogeneous equation, normalized such that 
R°°[r — > oo) ~ (ioor) 3 exp(iu)r*) . 

The most convenient way of integrating Eq. ( |1.1| ) is by means of a Green's function, 
which should be chosen so as to incorporate the specified boundary conditions. The standard 
theory p4] suggests that the appropriate solution is 



poo/ \ rr T} H (r\ r°° 

R( r ) = fLALl / p 2 ( r ')T(r')R H (r')dr' + / p 2 (r')T(r')R°°(r') dr', (1.2) 

Wt i2M Wt Jr 

where Wt is a constant, equal to the conserved Wronskian of R H (r) and R°°(r). The 
misconception is precisely that Eq. (|1.2|) must be the desired solution. In fact it is not: 
In order for the standard theory of Green's functions to guarantee that R(r) as given by 
Eq. (|1.2|) satisfies the specified boundary conditions, the integrals must converge. This is 
not the case here. Because R H (r) and R°°(r) both have a component growing as r 3 when 
r — > oo, and because p 2 (r)T(r) falls off only as r~ 3//2 , those integrals diverge when r — > oo. 
The claim that Eq. ( |1.2|) enforces the specified boundary conditions is therefore unjustified, 
and in fact, is wrong. 

The method adopted in this paper for regularizing the Teukolsky equation goes as follows. 
For simplicity, I consider the specific case of a particle falling radially into a Schwarzschild 
black hole, having been released from rest at infinity. Generalization to other situations 
should be straightforward. 

Instead of the ill-defined particular solution ( |1.2| ), the starting point of this analysis is 
the most general solution to Eq. ( |1.1|) , which is written as 
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B+ [ p 2 (r / )r(r , )i? x> (r / )dr / l, 

(1.3) 

where a, b, A, and B are constants. This represents the gravitational radiation generated 
by the infalling particle, plus the free radiation that was initially present in the spacetime. 
Boundary conditions must be imposed to eliminate this free component. 

Now, because the integrals are ill defined when r — > oo, these conditions cannot be 
imposed immediately. Instead, integrations by parts are carried out, which are specifically 
designed to make the integrals well behaved. The boundary terms at r' = a and r' = b are 
then absorbed into the constants A and B, while the boundary terms at r' = r are carefully 
kept. At this stage one is still dealing with the most general solution to the Teukolsky 
equation, and this is a perfectly valid starting point for the discussion of boundary conditions. 
Since the integrals are now well behaved, there is no obstacle in setting a = 2M and b = oo. 
The new constants A and B are then chosen so that the solution satisfies the specified 
boundary conditions. By this procedure, regularization of the Teukolsky equation involves 
no unjustified manipulations nor modifications to the equation, and is entirely natural. 

The idea described in this paper appears to be completely trivial, and it is indeed ironic 
that such a solution to the problem of divergent integrals did not come forth much sooner. 
However, the work required to carry out the procedure is not, in itself, entirely trivial. The 
rest of the paper is devoted to a detailed presentation. It is organized as follows: 

For the purpose of integrating by parts, it is convenient to write the Teukolsky functions 
R H (r) and R°°(r) in terms of the related, and better behaved, Regge- Wheeler functions 
X H (r) and X°°(r). These functions, and the transformations relating them, are described 
in detail in Sec. II, which establishes many results used later on. 

Regularization of the inhomogeneous Teukolsky equation is carried out in Sec. Ill for the 
specific case of a particle falling radially into a Schwarzschild black hole. Generalization to 
other cases should proceed along similar lines, but I shall not pursue this here. 

The main results of this paper are summarized in Sec. IV. 

Notation. The following symbols appear frequently throughout the paper: / = 1 — 2M/r, 
r* = r + 2Mln(r/2M - 1), d/dr* = fd/dr, Q = iMw, z = (iwr) -1 , and x = (r/2M) 1 / 2 . 



R(r) 



R~{r) 



A+ f p 2 {r')T{r')R H {r')dr' 

J a 



+ 



K 



II. REGGE- WHEELER, TEUKOLSKY, AND CHANDRASEKHAR 

This section is devoted to the derivation of various results which will be used in the 
following section. Specifically, Sec. II A contains a discussion of the Regge- Wheeler equation 
33|, and a study of the asymptotic behavior of its solutions near r = 2M and near r = oo. 



Section II B does the same for the homogeneous Teukolsky equation The Chandrasekhar 



transformation , which relates a solution to the Regge- Wheeler equation to a solution of 
the homogeneous Teukolsky equation, is the topic of Sec. II C. Finally, the functions X H (r), 
X°°(r), R H (r), R°°(r), and the relations between them, are introduced in Sec. II D. 
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A. Regge- Wheeler equation 



The Regge- Wheeler equation |33|] compactly describes a metric perturbation of frequency 
uj and spherical-harmonic indices £ and m. It reads 

-n 

>X(r)=0, (2.1) 



d 2 2 rW + 1) 6M 
+ u -f 



where f — 1 — 2M/r and d/dr* = fd/dr, so that 

r* = r + 2M ln(r/2M - 1). (2.2) 

For two linearly independent solutions -Xi(r) and X 2 (r), the conserved Wronskian is given 
by 

W RW (X 1 ,X 2 ) = f{X x X' 2 - X 2 X[), (2.3) 

where a prime denotes differentiation with respect to r. Because the Regge- Wheeler equation 
is real, if X±(r) is a solution, then X 2 (r) = Xi(r) is also a solution, linearly independent 
from the first; the overbar denotes complex conjugation. 

1. Asymptotic behavior near r = 2M 



It can be seen from Eq. ( |2.1|) that near r = 2M, the Regge- Wheeler function must behave 
as exp(— iur*) or its complex conjugate. To obtain more information, let 

X(r)=Y(f)e-^ r \ (2.4) 

Then a short calculation shows that the new function Y(f) must satisfy 

(1 - fff Y" + [(1 -/)(!- 3/) - 4Q) Y' - [£(£ + 1) - 3(1 - /)] Y = 0, (2.5) 



where Q = iMu and a prime denotes differentiation with respect to /. Equation (pO|) can 
be integrated by by writing 

oo 

Y(f) = 1 + E «n/ n , (2.6) 

n=l 

where the normalization was chosen arbitrarily. Substituting this into Eq. (|2.5|) and setting 
each term of the resulting series to zero gives 

£{£ + !) -3 



a-2 



1-4Q ' 
(£-!)£(£ +!)(£ + 2) - 12Q 



(2.7) 



4(1-2Q)(1-4Q) 

The higher-order coefficients can be obtained from the recurrence relation 

n(n - 4Q) a n = [2(n - l)n + £(£ + 1) - 3] a n _i - (n - 3)(n + 1) a n _ 2 , (2.8) 
which holds for n > 3. 
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2. Asymptotic behavior near r = oo 

Near r = oo, the Regge- Wheeler function must also behave as exp(— iuor*) or its complex 
conjugate. As before, let 

X(r) = Y{z)e~ lwr \ (2.9) 
where z = (iur)^ 1 . This leads to the differential equation 

z 2 (l - 2Qz) Y" + 2(1 + z - 3Qz 2 ) Y' - [£{£ + 1) - QQz] Y = 0, (2.10) 
where a prime denotes differentiation with respect to z. Equation ( |2.10|) is solved by writing 



Y(z) = l + J2Kz n , (2.11) 



71=1 



where once again the normalization was chosen arbitrarily. This gives 

b x = \£{£ + l), 



(2.12) 



b 2 = l[(£-l)£{£ + l)(£ + 2)-12Q}. 
All other coefficients can be obtained from the recurrence relation 

2nb n = (£ - n + 1)(£ + n) 6 n _i + 2(n - 3)(n + 1)Q 6 n _ 2 , (2.13) 
which holds for n > 3. 

B. Homogeneous Teukolsky equation 



The homogeneous Teukolsky equation 23 compactly describes a curvature perturbation 



of frequency uo and spherical-harmonic indices i and m, and is given by 

R(r) = 0, (2.14) 



-rP{ r )-r +p 2 (r)U(r) 
ar ar 



where 



(2.15) 



C/( r ) = f- l [(urf -4iu(r-3M)] - (£-l)(£ + 2). 
If i?i(r) and i?2( r ) are two linearly independent solutions, then 

W T (i*i, Ra) = P (RiR' 2 - R2K), (2.16) 

where a prime denotes differentiation with respect to r, is the conserved Wronskian. It 
should noted that the Teukolsky equation is complex; the complex conjugate of a solution is 
therefore not a solution. For this reason, an ingoing mode of the Teukolsky equation, which 
is proportional to exp(— iur*), must be distinguished from an outgoing mode, proportional 
to exp(iu;r*). 
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1. Ingoing mode: Asymptotic behavior near r = 2M 



It is easy to see from Eq. ( |2.14 ) that near the horizon, an ingoing mode must behave as 
f 2 exp(— iur*). This motivates the substitution 

R(r) = fS i *(f)e^ r \ (2.17) 

which implies 

= (1 - /) 3 / S' m " + (1 - f){7f 2 - 10/ + 3 - AQ) S' m ' 

-[8f -{f + £ + 14)/ + f + £ + 6 + AQ] S m . (2.18) 

Here, a prime denotes differentiation with respect to /. A solution can be obtained by 
writing 

oo 

S^(f) = l + J2Pnf n , (2-19) 



n=l 



which gives 



in _ £(£ + 1) + 6 + AQ 
Pi — 



P2 



3-4Q 

+ 1) + 4] [£{£ + 1) + 18] + 4[2£(£ + 1) + 33}Q 



(2.20) 



8(1-Q)(3-4Q) 
All other coefficients can be obtained from the recurrence relation 

n{n + 2 - AQ) $ = [n(3n + 4) + £(£ + 1) - 1 - A(n - 2)Q] 

- [n(3n + 2) + {£-!){£ + 2)] p™_ 2 + (n 2 - 1) p™_ 3 , (2.21) 

which is valid for n > 3 (with pJJ 1 = 1). 

2. Outgoing mode: Asymptotic behavior near r = 2M 



Equation (|2.14|) implies that near the horizon, an outgoing mode must behave as 



exp(— iur*). This leads the substitution 

R(r) = S out (f)e iujr \ (2.22) 

which implies 

= (1 - /) 3 / S out " + (1 - /)(3/ 2 - 2/ - 1 + AQ) S out ' 

- [(£ -1)(£+ 2)(1 - /) + 12Q) S out = 0, (2.23) 

where a prime denotes differentiation with respect to /. This is integrated by writing 
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S out (/) = l + £K U 7 n , (2.24) 



n=l 



which gives 



7iOUt _Ji-m + 2) + l2Q 
Pi — 



P2 



out 



1-4Q 



(£ - + 1)(£ + 2) + 12[2£(£ + 1) - 3]Q + 192Q 2 



(2.25) 



8Q(1-4Q) 

All other coefficients can be determined with the recurrence relation 

n(n - 2 + 4Q)p° n nt = [(n - l)(3n - 5) + (£ - l){£ + 2) + 4(ra - 2)Q] 

_ [(n _ 2 )(3n - 4) + (£ - l){£ + 2M U -2 + (n - 3)(n - 1) p£ s , (2.26) 

which is valid for n > 3 (with pg ut = 1). 



3. Ingoing mode: Asymptotic behavior near r = oo 

Near infinity, an ingoing mode of the Teukolsky equation must behave as r _1 exp(— iur*), 
which suggests the substitution 

R(r) = zS' m {z)e- iujr \ (2.27) 

where z = (iur)' 1 . The homogeneous Teukolsky equation then implies 

= (1 - 2Qzfz 2 S' m " + 2(1 - 2Qz)(l + 3z- 5Qz 2 ) S' ml 

+ [-{£ - 2){£ + 3) + 8Q + 2{£ 2 + £- 9)Qz + 12Q 2 z 2 ] S' m = 0, (2.28) 

where a prime denotes differentiation with respect to z. Once again the solution is written 
as a series, 

oo 

5 in (^) = l + EO n , (2-29) 

n=l 

and substitution yields 

q ? = l[(£-2)(£ + 3)-8Q], 

(2.30) 

Q2 = |{(^ " 3)(^ - 2)(£ + 3)(£ + 4) - A[A£(£ + 1) - 39]Q + 32Q 2 }. 

The recurrence relation 

2nC = [-n(n + 3) + (£ - 1)(£ + 2) + 4(n - 3)Q] C_ : 

+ 2[n(2n + 1) - £{£ + 1) - 1]Q <£_ 2 - 4n(n + 2)Q 2 C_ 3 , (2.31) 

valid for n > 3 (with = 1), gives the remaining coefficients. 
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4- Outgoing mode: Asymptotic behavior near r = oo 



Finally, an outgoing mode of the Teukolsky equation behaves as r 3 exp(iur*) near infinity. 
This leads to 

R(r) = z- 3 S out (z)e iujr *. (2.32) 

Equation (|2.14j) then implies 

(1 - 2Qz)z 2 S out " -(2 + 2z- 6Qz 2 ) S ouU - [{£ - 1){£ + 2) + 6Qz] S out = 0, (2.33) 
where a prime denotes differentiation with respect to z. This is integrated by substituting 

oo 

S out ( z ) = 1 + Y,Vn Ut z n (2-34) 

n=l 



into the differential equation, which gives 

?r = -i(i -ix*+ 2), 



(2.35) 



g° ut = \[{l + W + 2) -12Q]. 

The other coefficients are generated by the recurrence relation 

2n qT = [in - 4)(n - 1) - (I - 1)(£ + 2)] - 2(n - 3)(n - 5)Q g^, (2-36) 
which is valid for n > 3. 

C. Chandrasekhar transformation 

In 1975, Chandrasekhar [35]| proved the following theorem: If X(r) is a solution to 



the Regge- Wheeler equation ( |2.1| ), then there exists a linear differential operator C such 
that R(r) = CX(r) is a solution to the homogeneous Teukolsky equation (|2.14j) . The 
Chandrasekhar transformation is given by C oc r 2 / '£/ _1 £r, where £ = fd/dr + zu;. Since 
X(r) satisfies a second-order differential equation, C can also be written in first-order form 

as 



C= (iur) {2(1 -3M/r + iojr)r£ + /[£(£+ 1) -6M/r]}. (2.37) 

The constant of proportionality was chosen arbitrarily. 

Using the results derived in Sec. II A, simple manipulations are required to prove the 
following statements. First, concerning asymptotic relations near r = 2M: 

If X ~ e~^ r * , then CX ~ 4Qa 2 / 2 e" iLjr * . (2.38) 

If X ~ e^ r *, then CX 8Q 2 (1 - 4Q) e iwr *. (2.39) 

Next, concerning asymptotic relations near r = 00: 

If X ~ e~ lwr \ then CX ~ 2b 2 (zwr) - ^-***. (2.40) 
If X ~ e iwr * , then CX ~ 4 (iwr) 3 e^ r * . (2.41) 

The constants a 2 and 62 are given by Eqs. (|2.7|) and ( |2.12| ), respectively. As they must, 
the asymptotic relations found here for CX(r) agree with the relations derived for R(r) in 
subsection B. 
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D. Linearly independent solutions 

1. Asymptotic relations 

Of all the solutions to the Regge- Wheeler equation, two are preferred. The first describes 
gravitational waves which are purely ingoing at the black-hole horizon, and is denoted X H (r). 
The other describes waves which are purely outgoing at infinity, and is denoted X°°(r). 
These solutions satisfy the asymptotic relations 

where A m and *4 out are constants, and 

{€ tuJr T — > OO 
^gin g-iwr* _|_ gout gi^r* r y 2j\/f ' (2.43) 

where B m and B out are also constants. These solutions are linearly independent. Evaluation 
of their Wronskian in the limit r — > oo indeed reveals that 

W RG (X H , X°°) = 2iuA m , (2.44) 

where H^rw was defined in Eq. ( j2.3|) . 

Acting with the Chandrasekhar transformation, 

R H (r) = X H CX H (r), 

(2.45) 

R°°(r) = X 00 CX°°(r), 

where \ H and x°° are normalization constants, returns solutions to the homogeneous Teukol- 
sky equation possessing the same physical interpretation. These are normalized so that 

F>h ( \ / f 2 e~ iujr * r^2M 

^ ~ \ Q [n {icur)- 1 e- i " r * + Q out (iujrfe iuir * r -> oo ' ^ ' 

where Q m and Q out are constants, and 

noo W / (^r)V" r * r -> oo 

JT V J ~ S -pin y2g-iojr* _|_ pout gja;r* ^ ^ 2ikf ' V '/ 



where P m and P out are also constants. It follows from Eqs. ( |2.38|) and fl2.41| ) that this 
normalization is obtained by choosing 

h __ (1-2Q)(1-4Q) 

4Qa 2 Q[(£-l)^ + l)(£ + 2)-12Q]' 

(2.48) 

X 4 

These solutions are also linearly independent, and 

W T (R H , R°°) = -2^ 3 Q in , (2.49) 
where Wt was defined in Eq. ( [2.161) . 
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2. Relations among constants 



The constants .4 in < out , B' m ' ou \ Q in > out , and p in > out are not all independent. Various rela- 
tions among them are easily derived. 

In Eq. ( |2.44 ), the Wronskian Wrw(X , X°°) was evaluated near r — > oo. It can also be 
evaluated near r = 2M. Since the two values must agree, we have 

B out = A m . (2.50) 

Similarly, constancy of Whw(X h ,X°°) implies 

while constancy of W^sn{X h , X H ) gives 

|^in|2 _ |^out|2 = ^ ^ 

which expresses global conservation of energy. 

Additional relations are a consequence of the Chandrasekhar transformation. Combining 
Eqs. (ggg)-(ggg), ( ggg ), p^6| ), and ( ggg ) reveals that 

in _ J2_ , in (l-2g)(l-4g) in 
W 2Qa, AQ 

gout _ ^ ou * 

Qa 2 

Similarly, combining Eqs. Q), ( ggg) , (^43|) , ( ^45|) , ( ^47|) , and ( ^48|) gives 

7> in = Qa 2 £ in , 

p out = -2Q 2 (1 -AQ) B out . 
Finally, combining Eqs. (^50|) , flT5lD , ( ^53|) , and ( ^54|) yields 

pout _ Qin 

1 - 2Q ' 

V in = Q 2 \a 2 \ 2 Q out . 
This last equation does not follow easily from Wronskian relations. 



(2.53) 



(2.54) 



^2.55) 



III. REGULARIZATION OF THE TEUKOLSKY EQUATION 

I now proceed with the regularization of the Teukolsky equation. The source function 
T(r) is constructed in Sec. Ill A for the specific case of a particle of mass [i released from 
rest at infinity and falling with zero angular momentum into a Schwarzschild black hole of 
mass M. The general solution to the inhomogeneous Teukolsky equation is also displayed 
here. In Sec. Ill B the regularization procedure is carried out. Then the behavior of the 
regularized solution is examined near r = 2M in Sec. Ill C, and near r = oo in Sec. Ill D. 
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A. Source term and general solution 



The source function is easily constructed by following the steps spelled out in Poisson 



and Sasaki |23]. When the motion is purely radial, it is given by 

T(r) = 2y/(£ -!)£(£ + !)(£ + 2) r 4 f dt'dQ' T aP n a n p Y tm (6' , 0') e iuJt ' , (3.1) 



where dQ' = d cos 9 



I AMI n a = i, 



1, /, 0, 0) is a null vector pointing outward, Y^ m are the usual 



spherical harmonics, and T af3 is the particle's energy-momentum tensor, 

T a/3 (x') =/i [ dru a u^5[x' -x(r)]. 



(3.2) 



Here, x' represents an event in spacetime, labeled by the Schwarzschild coordinates 
(t',r',8',(f)'), and x(r) represents the particle's world line, with four-velocity u a = dx a /dr, 
where r is proper time. In Eq. (|3.2| ), the ^-function is normalized so that / S(x) v / — g d 4 x = 1, 
where g is the determinant of the metric. 

The geodesic equations for radial motion reduce to 9 = <fi = and 



dt 
dr 



7 V2M 



which integrates to 



t(r) 



-2M[ |x 3 + 2x + ln 



x — 1 



X 



(3.3) 

(3.4) 
1/f and 



where x = (r/2M) 1 / 2 . The four- velocity has non- vanishing components it* 
u r = —1/x. 

To obtain the source, Eq. (|3.2|) is first integrated with respect to dr, which returns the 
factor [iu a u^ jr 2 u r multiplying 5[t' — t(r)]6 (cos & — l)5((f>'). Contractions with n a are then 
taken and the result is substituted into Eq. (|3.1|) . After simplification, the result is 



where 



G 



p 2 (r)T(r) = Gg(r)e iu)t{r \ 



12 r(£-i)£(e + i)(£ + 2)(2e+i)] 1 / 2 



8M 2 



47T 



(3.5) 



(3.6) 



if m = 0, and G = otherwise (m is the spherical-harmonic index; that only modes with 
m = contribute to the full perturbation reflects the axial symmetry of the problem). Also, 



9{r) 



x(x + l) 2 



(3.7) 



The general solution to the inhomogeneous Teukolsky equation is obtained by substitut- 
ing Eq. (|3.5| ) into (|1.3|). The result is 
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+ R H (r 



J a 

B+ [ b £(r / ).R °(r')e iwt(r ' ) dr' 



(3.8) 



Here, W T = W T (R H ,R°°) = -2iu 3 Q} n , and a, b, A, and B are constants. 

Our task now is to see to it that the boundary conditions — waves ingoing at the horizon 
and outgoing at infinity — are properly imposed. In fact, there is no difficulty in demanding 
the correct behavior at the black-hole horizon. A short calculation indeed reveals that when 
a = 2M and A = 0, the first term of Eq. ( ft.8|) is 0(f 3 ). Since the second integral is finite, 
this ensures that R(r) oc R H (r) when r — > 2M, as required. Unfortunately, the behavior at 
infinity cannot so easily be controlled. This is because both integrals diverge when r — > oo, 
due to the fact that g(r) = 0(r" 3 / 2 ) while R H >°°(r) = 0(r 3 ). Clearly, the solution Q) 
must be regularized before an attempt is made to impose the correct boundary condition at 
infinity. 

B. Regularization 

I begin by defining the integrals 

2) = ^ g(r)R A (r)e iujt{r) dr, (3.9) 

where the index A stands for either U H" or "oo". For the purpose of regularization, R A (r) 
is conveniently expressed as 

R A (r) = X A CX A (r), (3.10) 
where C is given in Eq. ( |2.37| ). Substitution gives 

J A (l,2) = X A (/conv + /div), (3.11) 



where 



/conv = J" r conv (r)e^X» dr, (3.12) 
Idiv = [ T div (r)e l ^CX A (r) dr, (3.13) 



and 



r con v(r) = {iur)f[£{£+l) - 6M/r}g(r), (3.14) 

T div (r) = 2{iur){l - 3M/r + iur)rg(r). (3.15) 

As the names indicate, J conv is convergent when r — > oo, since r conv (r) = 0(r -1 / 2 ), while 
Jdiv is divergent, since T div (r) = 0(r 3//2 ). 
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Regularization of Idw can be achieved by integration by parts. To identify what must be 
done, consider the alternative form 



Idiv = j* T div e iuJt CX A + — (he iuJt CX J 
= I' div + boundary terms, 



dr - he iwt CX A 



(3.16) 



where h(r) is a function to be determined. After simplification, the new integral becomes 



4v = [ ^\r Ahl cx A + r; onv ^) dr, 



where 



r / dh 

T ** = Tr +lU3 



dt 1' ,. 



and 



r' 



+ 1) - 6M/r] — . 



(3.17) 



(3.18) 



(3.19) 



The function h(r) must be chosen so that I' div is well behaved when r — ► oo. 

The divergence of Idiv is caused by the bad behavior of r d i v (r). Happily, its contribution 
to J div can be removed by simply setting T' div = 0, which gives a differential equation for 
h(r). One solution to this equation is 



h(r) = 8M 



2 1 + x + 2Qx 3 



1 + x 



(3.20) 



where x = (r/2M) 1 / 2 . Substituting this into Eq. (|3.19|) then reveals that r^. onv (r) = 0(r _1 ), 
which implies that I div is indeed well behaved. 

Regularization has thus been achieved. Combining Eqs. (|3.11|) , ( p,12j ), (|3.16|) , ( p,17|) , 
and (|3.19|) gives 



i A a,2) 



X A [ g{r)e lut{r) X A {r) dr - X A h{r)e lujt{r) CX A {r) 



where 



X 



g(r) 



J A (l, 2) + boundary terms 



2(1 + Qx z 



x' 1 



1) -6M/r). 



(3.21) 



(3.22) 



This result will now be put to use. 

With the notation introduced above, Eq. ( j3.8|) reads 

G 



R( r ) = __{ J R-( r )[A + / H ( G ,r)] + J R H (r)[ J B + / 00 (r,6)]}. 



(3.23) 
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The regularized version of this equation is obtained by substituting Eq. ( |3.21|) . The boundary 
terms at r' = a and r' = b can be absorbed into the constants A and B by making the 
replacements 



A -> A + x H h{a)e luJt(a) CX H {a), 

B ->B- x^hiby^CX^ib). 
The boundary terms at r' — r combine to give 

J(r) = he^^B^CX 00 - x H R°°CX H ) 



(3.24) 



(3.25) 



This can be simplified by expressing R H '°°(r) in terms of X H >°°(r), as in Eq. ( p,10|) . After 
simplification, Eq. ( |3.25| ) becomes 



J( r ) = x H X°°W RW (iujr)f[£(£+ 1) - 6M/r]h(r)e 



iujt(r) 



(3.26) 



where W RW = W RW (X H ,X°°) = 2iuA m . 
Finally, gathering the results yields 



R{r) 



G_ 
Wt 



R°°(r) 



A + X H J H (a } r) 



+ R H (r) 



B + x °°J°°(r,b) +J(r) 



(3.27) 



where J(r) is given by Eq. (ET26ft , and 

^(1,2)= /% 



(r')e^ (r,) X A (r') d/, 



(3.28) 



as was first written in Eq. ( |3.21| ). 

The function R(r) is now expressed in terms of integrals that are well behaved when 
r — > cxd. There is therefore no obstacle in setting 



a = 2M, 



oo, 



(3.29) 



which I shall do from now on. It can then be verified that with this choice, the replacements 
of Eq. (|3.24j) take the form A — > A + 0, and B — ► B + oo. The infinite shift in B reflects 
the fact that the original expression for R(r), given by Eq. ( |1 . 2p , was not well defined. This 
shows the importance of the procedure carried out here: the integrals must be regularized 
before b is set equal to infinity. The fact that B is then shifted by an infinite amount is 
of no consequence: Since Eq. ( |3.27| ) is a general solution to the inhomogeneous Teukolsky 
equation, as can be verified by direct substitution, this equation is a perfectly valid starting 
point for the discussion of boundary conditions, to which I turn next. One might just as 
well forget how Eq. ( p. 27] ) was derived, and proceed afresh from here. 
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C. Behavior near r = 2M 



The behavior of R(r), as expressed by Eq. ( p.27| ), must now be examined near r = 2M, 
to ensure that it correctly represents purely ingoing waves at the black-hole horizon. The 
constants A and B must therefore be chosen so that R(r — > oo) ~ (constant)/ 2 exp(—iujr*) . 

Our first task is to evaluate J H (2M,r) in the limit r — > 2M. Because R°°(r) = 0(f°), 
this calculation must be carried out to second order in /. On the other hand, only the 
leading-order term in J°°(r, oo) is required for the calculation; this is simply given by 
J°°(2M, oo). Finally, J(r) will have to be computed, also to second order in /. 



1. Evaluation of J H (2M,r) 
The results of Sec. II A imply that near r = 2M, X H (r) can be written as 

X H (r) = Y(f)e- iujr \ (3.30) 



where Y(f) = 1 + a\f + a 2 f 2 + 0(f 3 ). The coefficients a\ and a 2 are given by Eqs. 
Substituting this into Eq. ( |3.28| ) returns an exponential factor of the form exp(iuu), where 
u(r) is defined by 

u(r) = t(r) - r* = -4M[|x 3 + \x 2 + x + ln(x - 1)] , (3.31) 
where x = (r/2M) 1 / 2 . Changing the integration variable, what must be evaluated is 

rH,^,, f°° 9f Y „w 



J H (2M,r)=J j^e™ u du', (3.32) 

where the integrand is considered to be a function of u' . 

To compute the integral, Eq. ( |3.31| ) must first be inverted in order to express x as a 
function of u. While this cannot be done exactly in closed form, what is required here is a 
result accurate only to second order in / = 0(x — 1) = 0{e~ u ^ M ). Equation ( |3.31| ) implies 

U = ex p("^ " y) = (* ~ 1) + 3(* - l) 2 + 0[(x - l) 3 ], (3.33) 
which can be inverted to give 

x- 1 = U -3U 2 + 0(U 3 ). (3.34) 

Next, the integrand is expanded in powers of x — 1, and Eq. (|3.34j) is used to write this in 
terms of U. Integration is then straightforward. The final result must be expressed as an 
expansion in powers of /. For this purpose one uses U = \f + |/ 2 + 0(/ 3 ), which follows 
from x - 1 = \f + |/ 2 + 0(/ 3 ). The final result is 

J H {2M, r) = 4M^/ + ii 2 f 2 + 0(f)]e^\ (3.35) 

where 
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(l + Q)[£(£ + l)-3] 

(3.36) 

2(l 4 + 2£ 3 - U 2 - 61 + 12) + (2l 4 + 4l 3 + 111 2 + 91 - 63)Q + [61(1 + 1) - 42]Q 2 

4(1-2Q)(1-4Q) ' 

The evaluation of J H (2M, r) is now completed. 

2. Evaluation of J(r) 



This calculation is quite straightforward. Equation fl3.26| ), with h(r) given by Eq. ( |3.20| ) 



can immediately be expanded in powers of /. The result is 

J(r) = 16M 2 Qx H X°°WW[vi/ + + 0(/ 3 )]e^, (3.37) 



where 



u 1 = (l + Q)[£{e+l)-3], 



v 2 = £(£+l) + ~[3£(£+l)-5]Q. 



(3.38) 



5. Evaluation of R(r) 

Expansions in powers of / have been obtained for J H (2M,r) and J{r). These are 
supplemented by the results of Sec. II B and D, which imply 

i?°°(r) =V in [f 2 + 0(f)]e~ iujr ' +V out [l+p° 1 ut f + 0(f)]e iujr \ (3.39) 

where p° nt is given by Eq. ( |2.25| ). We also have R H (r) = [f 2 + 0(f 3 )} exp(— iur*) and 
J°°(r, oo) = J°°(2M, oo) + 0(f). 

Substituting all this into Eq. (|3.27| ) gives R{r) as an expansion in powers of /, which 



contains terms of order /°, /, as well as the allowed terms of order f 2 and higher. Each 
term will be discussed in turn. 

The term of order f° can be eliminated by setting 

A = 0, (3.40) 

which will be done from here on. 

The term of order / can be simplified using Eqs. ( gig ), (^48|) , (|2753|) , (|2~55| ), ( ^36|) , 
and ( 3.38| ). As it must, it vanishes identically. 



The term of order f 2 survives, and contains two contributions. The first is proportional 
to exp(icjt) and comes from J H (2M,r) and J(r); the other is proportional to exp(— iur*) 
and comes from R H {r). The first contribution is simplified using the same equations as 
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before, in addition to Eq. (gjg ); the result is AMQ(1 + Q)A in f 2 e iuJt . This is then combined 
with the second contribution by observing that t = —r* — 2M(| — 2 In 2) + 0(f). 



After simplification, the final result is that near r = 2M, 
G 



fi + C + x 00 ^ °(2M,oo) f 2 e~ lujr * + 0(f) 



(3.41) 



as required. It is recalled that J°°(2M, oo) was defined in Eq. (|3.28|) , and the constant C is 
given explicitly by 



C = AMQ(1 + Q)A m exp 
The constant B will shortly be set to zero. 



-2Q(--21n2 



(3.42) 



D. Behavior near r = oo 



It is much easier to extract the behavior of Eq. (|3.27| ) near r = oo. I begin with the 
computation of J°°(r, oo). In this limit, Eq. ( |3.22| ) reduces to 



9{r) 



2Q£(£ + I) 



x 



l + 0(x- 



while Eqs. (|2]|) and ( [DID imply 

X°°(r) = [1 + 0(ar 2 )] e^ r *. 
After substitution into Eq. (|3.28|) , and a change of integration variable to 



v (r) = t(r) + r* = -4M[|x 3 - \x 2 + x - ln(x + 1)], 



(3.43) 



(3.44) 



(3.45) 



one arrives at 



J°°(r, oo) = 2£(£+l)Q / — 1 + 0(x^ 



v 1 



oo X 



e l0JV dv' 



+ l)Q 2 z[l + 0(z 1 l 2 )}e l > 



(3.46) 



where z = (iuir)~ l . 

The computation of J(r) is also straightforward. Equations (|3.20|) and ( p.26|) immedi- 
ately give 



J(r) = 8M 2 /fMHl)^ l + 0(z 



1 



e i0Jt . 



(3.47) 



To finish the job, Eqs. (|3.46|) and (|3.47|) , together with the relations J H (2M,r) = 
J H (2M, oo) + 0(z 1 / 2 ) and R°°(r) = [z~ 3 + 0(z~ 2 )] exp(icur*), are substituted into Eq. ( &2~7\) . 
This gives 



R(r) = [ X H J H (2M, oo)[l + 0(z 1/2 )}R°°(r) + BR H (r) + other terms}, (3.48) 



18 



where the "other terms" are all proportional to exp[iu)t(r)], and are 0(z 2 ) or higher, and 
therefore much smaller than the dominant terms of order z~ 3 . Furthermore, because t(r) ~ 
-(4M/3)(r/2M) 3 / 2 when r — > oo, their phase increases much more rapidly than r* ~ r, 
which means that the "other terms" cannot be combined into a term proportional to R H (r), 
which would represent a free, initially incoming, gravitational wave. The only such term 
present in R(r) is BR H (r), and the requirement that waves must be purely outgoing at 
infinity dictates 

B = 0. (3.49) 

The final result is that near r = oo, 

R(r) = — X H J H (2M, oo) (iuurfe luJT * + 0(r 5/2 ), (3.50) 
as required; J H (2M, oo) was defined in Eq. ( |3.28| ). 



IV. SUMMARY AND CONCLUSION 



I now summarize. The regularization of the Teukolsky equation was successful. The 
procedure consisted of two stages. In the first stage, the most general solution to the 
inhomogeneous Teukolsky equation was written in terms of integrals that are well behaved 
when r — > oo. This was accomplished in Eq. (|3.27|) . In the second stage, an ingoing- wave 
boundary condition was imposed at the black-hole horizon, and an outgoing-wave boundary 
condition was imposed at infinity. While the correct behavior at infinity could not be verified 
with the original form of the solution given by Eq. (|3.8|), it was quite straightforward to do 
so with the regularized form (|3.27|) . 

The regularized solution is written as 



where 



R(r) 



G 



X H J H (r)R°°(r) + X °° J°°(r)R H (r) + J(r) 



G 



8M 2 
W T = -2iu 



X 



X 



H 



3 V 1-2Q)(1-4Q) 
4Q 

(1-2Q)(1-4Q) 



A 11 



Q[(£ -!)£(£ + l)(£ + 2) - 12QY 
1 

4' 

g(r')e itJt( - r,) X H (r')dr', 



J H {r)- 

Jim 

/oo 
^(r')e^ (r ' ) X 00 (r') dr', 

J(r) = -2LU 2 x H x 0O A m rf[£(£+l) - 6M/r]h(r)e iuJt{r) . 



(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 
(4.8) 
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The functions X H (r) and X°°(r) are linearly independent solutions of the Regge- Wheeler 
equation ( |2.1|) , with normalizations determined by Eqs. ( 2.42|) and (|2.43|) . The constant A 111 
is also defined by these equations. Also, / = 1 — 2M/r, x = (r/2M) 1 / 2 , Q = iMu>, and 



9{r) 
h(r) 
t(r) 



2(1 + Qx 3 



x 1 



+ 1) -6M/r], 



8AP 



1 + x + 2Qx 3 
1 + x 



-2M( §ar 3 + 2x + ln 



(4.9) 
(4.10) 
(4.11) 



J 



Equation ( f4.1|) is obtained directly from (|3.27|) by setting a = 2M, b = oo, A = B = 0, 



J H (2M,r), and J c 



J°°(r, oo). Equation ( |4.2| ) is the same as ( |3.6| ). Equation 



( fl~3"D follows from ( |2.49| ) and fl2.53Q . Equations (0) and ( ^5|) are the same as ( p.48|) . 
Equations ( |4.6| ) and ( [4.7| ) are the same as ( |3.28| ). And finally, Eqs. ( |4.8| )-( [4TTT| ) are the same 
as (|3~26D , (1021) , (|3^0|) , and Q, respectively. 

Equation ( [4.1[ ) implies that near r = 2M, i?(r) behaves as 



i?(r) 



G_ 



2 2+4Q e- 10Q/3 MQ(l + g)^ in + x°° J°°(2M)1 / 2 e 



This follows from Eqs. ( |3.41|) and ( |3.42| ). It also implies that near r 

G 



(X), 



R(r) 



X H J H {oo) {luorfe 



(4.12) 



(4.13) 



which is the same statement as in Eq. (|3.50j ). This expression agrees precisely with the one 
derived by Simone, Poisson, and Will 0], who obtained it by "throwing away the infinite 
boundary term" . 

The radial function R(r), found here to be a solution of the inhomogeneous Teukolsky 
equation, represents a gravitational perturbation of frequency to and spherical-harmonic 
indices I and m. A better notation for it (which I did not adopt in order to keep all symbols 
simple) would be Re m (uj;r). The full perturbation is obtained by summing over all these 
modes. More precisely, the perturbation in the Riemann tensor caused by the infalling 



particle is represented by the complex function \1>4 |36| given by 



r J im 



—iu)t 



du. 



Here, x is the event in spacetime labeled by the Schwarzschild coordinates (t, r, 8, 
-2Yim(9 , 4>) are spherical harmonics of spin-weight —2 |37 



the functions 



(4.14) 

»), and 
For an axially 



symmetric problem, such as the one considered in this paper, modes with m ^ vanish 
identically, so the sum over m reduces to the single term m — 0. From ^(x) one may 
obtain many relevant quantities, such as the gravitational- wave field h^(x) and the fluxes 
of energy at the black-hole horizon and at infinity. Additional details are provided by 
Refs. 



23.3C 



The considerations of this paper were limited to the simplest case of an infall into a 
black hole: the hole was assumed to be nonrotating, and the particle was assumed to have 



20 



zero angular momentum and a vanishing initial velocity. There is, however, no reason to 
believe that the methods used here could not be extended to more complicated situations. 
Of course, the amount of labor involved, already considerable here, would increase, but there 
is no issue of principle. 

To conclude, I would like to stress the main message of this paper. The standard choice 
of Green's function for solving the inhomogeneous Teukolsky equation leads to divergent 
integrals, and contrary to naive expectations, fails to enforce the correct boundary conditions 
at the black-hole horizon and at infinity. The regularization procedure amounts to nothing 
more — and nothing less — than finding an adequate Green's function. Contrary to what 
may have been believed, there is nothing intrinsically wrong with the Teukolsky equation 
when dealing with non-compact source terms. 
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